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Abstract 

The existence of static, self-gravitating elastic bodies in the non-linear theory of elasticity is estab- 
^v^j . lished. Equilibrium configurations of self-gravitating elastic bodies close to the reference configuration 

have been constructed in [6] using the implicit function theorem. In contrast, the steady states considered 
| in this article correspond to deformations of the relaxed state with no size restriction and are obtained 

Q\ • as minimizers of the energy functional of the elastic body. 
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1 Introduction 



Steady states of bodies subject to their own self-induced gravitational field are at the center of interest 
in theoretical astrophysics. They describe different physical systems according to the specific type of 
material making up the body. The most popular examples are steady states of self-gravitating fluids, 
which describe equilibrium configurations of stars [10], and of kinetic collisionless matter (Vlasov matter), 
which are models for steady galaxies [7j. For these important examples there is a vast mathematical 
literature concerning the existence theory of static solutions, see for instance [SJ 1151 1161 1231 1241 1251 126| . 
In contrast, the existence of static self-gravitating elastic bodies has been scarcely investigated so far, 
despite the prominent role that these models play in the context of neutron stars physics [9]. In fact, 
the only result that we are aware of is [5] (and its extension to the general relativistic case Q]), where 
R. Beig and B. G. Schmidt prove the existence of static, self-gravitating elastic bodies near the reference 
configuration by using the implicit function theorem. In the present article we consider elastic bodies 
that correspond to deformations of the reference state without any a priori size restriction. The precise 
formulation of our results is given in Section [2] This Introduction continues with a general discussion on 
the problem of self-gravitating bodies in equilibrium. 

In the approximation in which thermal effects can be neglected, the mass density p(x) of isolated, 
self-gravitating matter in equilibrium satisfies the equation 

Divrj = pW, V(x) = - [ P ^ dy, i£l 3 , (1) 
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where a = u(x) is the Cauchy stress tensor and V — V(x) is the gravitational potential self-induced by 
the matter distribution. 

The system fT]) must be complemented by a constitutive law expressing the stress tensor in terms of 
the mass density. For simplicity we assume that this is given by an explicit equation a — a(x,p), with 
&ij(x, 0) = 0, for all x G R 3 and i, j = 1, 2, 3. For instance, for a barotropic fluid with pressure p = p(p), 
p(0) = 0, the stress tensor takes the form ay = —pSij, where <5y = 1 if j = i and 5y = otherwise. 

Definition (Single Body). The matter distribution is called a single body if there exists a non-empty 
open bounded connected set Q C R' ! with boundary d£l of zero Lebesgue measure such that p > 0, for 
x G Q and p = 0, for Q c . 

Note that our definition of single body does not require regularity of the boundary dtt. Moreover, 
since the mass density of a single body need not vanish on dQ, then p is in general only a weak solution 
of Q. We recall that a function p : R 3 ->■ [0, oo) is a weak solution of Q if p G L 1 nL 3/2 (R 3 ), a G L 1 (R 3 ) 
and the following integral equation holds: 

/ Tr(a • V</> T ) dx = - / / ( x - y ) . fa) dx dy, (2a) 

Jr3 Jm.3 J R 3 \x - y\* 



for all : R 3 ->• R 3 , <j> G C,!(R 3 ), (2b) 

where A ■ B denotes the dot product of the tensors A and B (i.e., the contraction of the last index of A 
with the first index of B), while A T and Tr(^4) denote the transpose and the trace of the second order 
tensor A, respectively (our convention for the gradient is (V/)ij = d x jfi, for all functions / : R 3 — ► R 3 ). 

For single body solutions, we may restrict the integrals in (f2]) over the support fl, and obtain that p 
is a weak solution of {TJ if the the mass density in the interior of the body solves the integral equation 

/ Tr(cr • Vr? T ) dx = - / / P( x )P(y) ( x -y). n (x) dx dy, (3a) 



for all 77 : Q -> R 3 , 77 G r>n(R 3 ), (3b) 

where X>n(R 3 ) is the trace of Cc(R 3 ) on f2, i.e., the space 

Dn(R 3 ) = {? 7 6 C^n) : 77 = 0| n , for some : R 3 -> R 3 , G C^R 3 )}. (3c) 

Since 2?n(R 3 ) C C 1 (?l) C i? 1 ^), a subset of solutions are those that solve (J3aj) for all 77 G H 1 ^). In 
the latter case, and assuming that dQ, is smooth, p is a weak solution of the following Dirichlet-type 
boundary value problem 

Diva = p/ 1 P{V \ 2 ^ V , dy, x&n, (4a) 

Jn \x-y\ 2 \x-y\ 



p = pb : cr(x, pt) • 71 = 0, x G dfl, (4b) 

where n denotes the exterior unit normal vector field on dQ. We also remark that compactly supported 
solutions of ([T]) might exist only for special domains Q. For instance, isolated, self-gravitating fluid bodies 
are necessarily spherically symmetric [171 119| and therefore they are all supported on a ball. 

Problem (|3]) is posed in the Euler formulation of continuum mechanics, which is the natural one 
for fluids. However for self-gravitating elastic bodies the problem is more naturally formulated in the 
Lagrangian picture, since the stress tensor of an elastic body depends on the deformation of the body; 
see [12] for an introduction to 3-dimensional non-linear elasticity. 

Let B be an open, bounded, connected subset of R 3 with smooth boundary and ip : B — ¥ R 3 be an 
injective function on B. Denote 

F = W, Cof F = F~ T det F (the matrix of cofactors of F). 

It is assumed that 7/) preserves orientation, i.e., det Vt/>(X) > 0, for all X £ B. The set B is called reference 
configuration (or material manifold) and identifies the body in a completely relaxed state, before any 



2 



deformation takes place. Upon a deformation, the "particle" X of the body is moved to the new position 
x = ip(X). If we denote by p re f : B — > (0, oo) the mass density of the body in the reference configuration, 
and by p : ip(B) — t (0, oo) the mass density in the deformed state, the conservation of mass entails 

The (first) Piola-Kirchhoff stress tensor is defined by 

T,(X)=a(X)-CoiF(4*(X)). (6) 
The change of variable x = ip(X), together with ([5]) and ((6]l turns equation (|4a|) into the form 

DivS = Pret f s p re i(Y) jfff,, *r, XSB. (7) 

By inverting ([5} and ((6]l , sufficiently smooth solutions of ((7)) are transformed into solutions of (|4afl in the 
domain Q, = ip(B). If the deformation ip is injective on the closure of B, then Q is a domain with smooth 
boundary. However, if injectivity is lost at the boundary, as it is often the case, there is no guarantee 
that the set fl possesses a smooth boundary, which is the main reason for introducing the definition of 
single body as above. We remark that interior injectivity prevents the interpenetration of matter, while 
injectivity up to the boundary prevents that even self-contact of the boundary of the body occurred. 

The body is said to be elastic if £ depends on the configuration ip only through the deformation 
gradient F = Vi/>. Denoting by M 3 the linear space of 3 x 3 real matrices, and by the subset 
thereof consisting of matrices with positive determinant, we then assume the existence of a function 
E : Bx M 3 — > M 3 such that E(X) = t(X,F(X)). For elastic bodies we may introduce a variational 
formulation of (Q as follows. Let a function w : B x — > [0, oo) be given and let us define the energy 
functional of the elastic body as 

PlA X M X, W) d X - \ f g I ^ff^ff, «* dY - (8) 

The function w is the stored energy of the elastic body. The first term in the energy functional is the total 
strain energy, while the second term is the gravitational potential energy of the body. A simple formal 
calculation shows that © corresponds to the Euler-Lagrange equation associated to the functional (|8j) 
for a Piola-Kirchhoff stress tensor given by 

* = (9) 

This suggests that, under appropriate assumptions on the growth and the regularity of w, one can obtain 
solutions of (O by proving the existence of minimizers of the functional / over a suitable space. The 
choice of the minimizing space is one of the main issues of the problem. We shall consider two different 
minimizing spaces, A (1) and A (2) , which in the pure elastic case were first considered in [3] [11]. The 
deformations in the space A^ 1 ' preserve the mass of any local subregion of the body, as well as the center 
of mass in the deformed state, but they are not in general injective functions, not even in the interior 
of the body. On the other hand, the deformations in the space A^ are homeomorphisms on B that 
preserve the shape of the body in the deformed state. 

The results concerning the existence of a minimum of the functional ([§]) are contained in Theorem 
[T] below. Once the existence of minimizers of the functional ([8]) has been proved, it is natural to ask 
whether they are solutions of (0. We are not able to give a positive answer to this question; even in 
the absence of self-gravitating interaction this is a major open problem in elasticity theory [5]. However, 
following the ideas introduced by J. Ball in [4], we are able to show that 

1. The spatial density associated to the minimizer tp^ 1 ' in the space A' 1 ' is a solution of ((3} in the 
region n' 1 ' = tp (B), and thus describes the interior mass density of an isolated, self-gravitating 
body in equilibrium; the boundary of fi^ 1 ' has zero Lebesgue measure, but we cannot say anything 
about its regularity. 

2. The spatial density p^ associated to the minimizer tp <,2 ' > in the space A^ solves equation (|4a[) in 
a distributional sense in the domain f2 (2) = i/j (2) {B); in this case, the boundary of Q.^ is regular, 
but since p' 2 ' need not satisfy ([2]), the body will not in general be isolated. 

The precise statement of our main results is given in Section [2] their proofs are to be found in Section [3] 
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2 Preliminaries and main results 

We assume that B is a regular domain, i.e., a non-empty, open, bounded, connected subset of K 3 with Lip- 
schitz continuous boundary; in particular, meas dB = 0. The mass density in the reference configuration 
p rc f : B —¥ (0, oo ) is assumed to satisfy 

Plc( eL°°{B), essinf p lel {X) = p > (10) 

and we denote 

M = ||Prcf|| L l. 

We write the energy functional of the elastic body as 

I{tp] = E str [V] + -Epot[V'], 

where E po t[tl>] is the potential energy, 

EpotW = -\ [ [ &i,{X,Y)p Ie t(X)p Iei (Y)dXdY, (11a) 

Z JB JB 

e *V> Y )= MX)-1>{Y)\ > (llb) 

while E s ti [V 1 ] is the total strain energy of the body, 

E stI [^]= f Prci (X)w(X,V^)dX. (12) 

JB 

Our first goal is to prove the existence of a minimizer to the functional I in two different spaces, which 
we denote A (1) and A (2) . We assume A^ C W 1 '^(B), for /3 > 3, and henceforth any function in the 
Sobolev space W 1 '^(B) will be identified with its representative in C°(B). Furthermore we assume that 
the elements of A^ satisfy det W > for almost all X G B. Since Vtp = a.e. on any set in which ip is 
constant (see for instance [131 Lemma 7.7]), we infer that \ip{X) — ijj(Y)\ > for almost all X, Y G B x B 
and therefore the function Q(X,Y) that appears in the potential energy (|1 la|) is well-defined almost 
everywhere in B x B. 

Before completing the definition of the spaces A^\ we need to list our assumptions on the stored 
energy function: 

(wl) Poly convexity: There exists a convex function w(X,-) : M 3 xM 3 x(0, +co) —¥ R such that w(X, F) = 
w(X, F, Cof F, det F), for all F £ M 3 and for almost all X G B; 

(w2) Regularity: The function w(-, F, H, 5) is measurable for all (F, H, 8) G M 3 x M 3 x (0, oo). 

In addition we impose a lower bound on the stored energy function which depends on the minimizing 
space considered. For F G M we denote |F| = s/Tt(F ■ F T ) (or any other equivalent matrix norm), 
while \a\ stands for the standard Euclidean norm when a is a real number. When working in the space 
AW 

we assume that the stored energy functions satisfies 
(w3) (1) There exist 

P>6, q>-^, s >^' ( 13 ) 
p — 1 p — 6 

a constant a > and a function h G ^(B) such that 

w(X,F) > a(|detfr s + |F| P + |Cof F\ q ) + h{X), 

for almost all X G B. 
When working in the space A^ we assume 
(w3) (2) There exist 

p > 3, q > 3, s > (14) 
q o 

a constant a > and a function h G i 1 (S) such that 

w(X,F) > a(|det^r s + \F\ P + | Cof F\ q ) + h(X), 

for almost all X G B. 
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Up to the choice of the exponents p, q, s, the assumptions listed above are standard in elasticity theory. 
For a throughout discussion on stored energy functions and examples that satisfy (or do not) the above 
conditions, we refer to [5] and |12l Chs. 3-4]. At the end of the present section we shall discuss briefly 
the class of Ogden materials. We also remark that there exist additional physical conditions that the 
stored energy functions must satisfy, which however are irrelevant for our analysis. For example, w must 
be normalized such that 

w(X,T) = 0, for all X £ B. (15) 
Next we complete the definition and prove some basic properties of the spaces A^- 1 ', A^ 2 K 

The space 

Let 

S^p — {x £ ip(B) : caxd{ip^ 1 (x)} > 1}, (16) 

where {ip~ 1 (x)} denotes the pre-image set of x £ ip(B) and cardt/ denotes the cardinality of the set U, 
and let a £ R 3 . Given a stored energy function that satisfies the properties (wl), (w2), (w3)^, we define 
the space as 

A w = {i>£ W hp (B) : £ 8t rM < oo, 

det X7ip > a.e. on B, 
meas = 0, 



p Iei (X)ip(X)dX = }. 

Since we can always find a constant matrix C £ M+ and a vector d £ R 3 (depending on a and p Te t ) such 
that ip(X) = C ■ X + d £ A {1) , then A (1) is not empty. 

A function ip that satisfies the condition meaaS^ — is often referred to as a.e. injective [T5]- The 
meaning of the last condition in the definition of A^ 1 ' will be clarified shortly. Before proving some basic 
properties of functions in the space A (1 \ let us note that 

S,p=ip{K 4 ,), where = {X £ B : 3 Y £ B, Y / X, such that tp(X) = ip(Y)} 

and we recall that the change of variables formula for Sobolev maps ip £ W 1,e (B), P > 3, is [2TJ Th. 2] 

f(ip(X)) detVip(X)dX = [ f(x)caxd{ip~ 1 (x)}dx, (17) 

for all measurable sets U C B and measurable functions / : ip(U) —¥ R. In particular, since functions in 
the space A^ 1 ' are a.e. injective we obtain 

f(ip(X)) det\7ip(X)dX = f f(x)dx, for all ip £ A^\ (18) 

Lemma 1. For all ip £ A^ the following holds: 

(i) ip maps null sets into null sets; 

(ii) ip is open; 

(hi) ip(B) is open, connected and satisfies 
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iP(B) = 4>{B), diP(B) C i>(dB); (19) 

(iv) meas ^ = meas dip(B) = 0; 

(v) The spatial density 

P * [X) - detVV^-H*))' 
is almost everywhere defined and positive on ip(B), and satisfies 

[ p tef (X)dX = [ p^(x)dx, for all ip £ A (1) . (21) 
In particular, ||p||ii(R3) = M, for all ip £ A^\ 
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(vi) The last condition in the definition of the space can be rewritten as 

/ xp^(x)dx = a, for all ip £ A w ; (22) 
Ji>{13) 

(vii) The potential energy can be expressed as 

E pot [tp] = epotDfy], (23a) 
where, denoting by lu the characteristic function of the set U, 

P(x)p(y) 



Pd, = PV^(B)» e P ot[p] = / 
Jut* 



k3 Je3 \ x - y\ 



■dxdy. (23b) 



Proof. The property (i) is valid in general for all functions tp £ W 1,fi (B), ft > 3, see [21]. To prove (ii), 
we first note that 

(|VV»r + (detV^)~ s ) dX < 



for all tp £ A"' , where p, s are given by (|13[1 . Let us introduce the outer distortion Ko of 

- detW 

By Holder's inequality, K £ L (B), P = ps/(p + 3s). Since /3 > 2, then (ii) follows by [H Cor. 1.10]. 
Next we prove (iii). By (ii), tp(B) is open; since continuous functions map connected sets into connected 
sets, it is also connected. The properties (|19p are proved in [121 Thm. 1.2-7] under the hypothesis that tp 
is injective on B, but the same result, with the same simple proof, holds under the (weaker) assumption 
that tp is open. As to (iv), applying the change of variable formula (|18l) we have 

det S7tp dX — dx — meas S,p = 0. 



Since det X7tp > a.e. in B, then meas K^,=0. Moreover, meas dtp(B) = follows by (i) and the inclusion 
dip{B) C ip(dB). To prove (v), let 

T il = {x£ tp(B) : det Vip{X) = 0, for some X £ {ip' 1 (x)}}. 

Since Ty, = tp({X £ B : det S/tp(X) = 0}), then measly = 0. As the spatial density ((20)) is well-defined 
and positive for x £ tp(B) \ (S^, UTy,), the first part of the claim follows. The identities (|2ip . (1221) and (|23|) 
are proved using (|18[) . □ 



We remark that by ([21[) . the deformations tp £ A^ preserve the mass of any subregion of the 
while by (|22[) , the deformations tp £ A^ 1 ' leave invariant the center of mass of the body in the deformed 
state. 

The space A^ 

For the definition of A {2 \ let a function £ : B ->• M 3 , £ £ W /1 ' P (B), p > 3, be given with the following 
properties: 

(zl) £ is one-to-one on 6, det V£ > a.e. in B and -E str [£] < oo; 
(z2) The boundary of the set ((B) is Lipschitz continuous. 

Note that under assumptions (zl)-(z2), C,(B) is a regular domain and that a sufficient condition for (z2) 
is that ( be a C 1 diffeomorphism on B. Moreover by [121 Th 1.2-8] we have 



C(B)=C(B), d((B) = adB). (24) 

Given a stored energy function that satisfies the properties (wl), (w2), (w3)' 2 ' and a function £ as above, 
we define the minimizing space A^ as 

A {2) = {ip£ W hp (B) : E stI [tp] < oo, 

det Vtp > a.e. on B, 
^(forie dB}. 

Since £ £ A^ 2 \ then A^ is not empty. 
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(r > 3); (25) 



Lemma 2. The properties listed in Lemma\j] are also valid for all tjj G 7l' 2 ' , and in addition the following 
holds: 

(A) V(B)=C(B); 

(B) if) is a homeomorphism of B onto ((B); 

(C) The inverse function tp~ belongs to W ' r (((B)), where 

r _ 9(1 + *) 
q + s 

(D) There holds 

^(S) = C(S), ^(B) = SC(B)- (26) 
/n particular, ip(B) is a regular domain. 

Proof. Since 

f (|cofV^| g + (detVtp)~ s )dX < oo, for ip G A (2) , 
Jb 

with q, s given by (|14p . Holder's inequality entails 

/ \(Vipy 1 \ r detVipdX = [ (detVip) 1 ~ r \CofVip\ r dX <oo, 
Jb Jb 

where r is given by (|25p . Thus the claims (A), (B) and (C) follow by [3J Th. 2]. Applying again |12l 
Th 1.2-8] we infer that the properties (|24p are also satisfied by ip, which gives (D). Finally, all the 
properties listed in Lemma [1] also hold for tp G A^ , because injective continuous functions are open 
(Invariance of Domain Theorem) and satisfy = 0. □ 

By (|26[) . it follows that the deformations in the space A^ preserve the shape of the body in the 
deformed state and that self-contact of the boundary does not occur. 

Main results 

Our first main result is the following: 

Theorem 1. Let p Ie { satisfy (|10[1 . Fix i = 1 or 2. Let w satisfy (wl), (w2), (wSf^ , and, for i = 2, let 
£ satisfy (zl), (z2). Then there exists G A^' such that 

= inf I[ip]. 

AW 

Moreover, the spatial density 

M 

P = fy>« 

satisfies p (i) G L 7 (^ W (B)), /or aZZ 1 < 7 < 1 + s. 

Next we study the relation between the minimizers of the energy functional and the mass density of 
static, self-gravitating bodies. To this purpose we need an additional assumption on the stored energy 
function: 

(w4) w(X, ■) is C 1 and there exists a constant K > such that 



dw 
8F 



<K(w(X,F) + l). 



Let us denote 

By Lemma [IJiii) and Lemma EJD), fi w is open, connected and meas af2 (i) = 0; ft (2) , in addition, 

has a Lipschitz continuous boundary and so is a regular domain. Define the Cauchy stress tensor 
<t w : fi w -» M 3 as 

o^Hx) - / ^(x,v^>(x)).v^W \ E _ , f ft» f27) 

a {X) -\ detV^(X) I ' L - prci 9F' (27) 
Our result concerning the equation solved by the minimizers is the following. 
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Theorem 2. In addition to the hypotheses of Theorem^ let w satisfy (w4)- Then a' 1 ' £ L 1 (r2^') and 
for i = 1 the following identity holds: 

Spot [Pa)] = / Tram da, (28) 

where 

- (l)ir — fflir 

= P IflM. °"(1) = °" If2(i)- 

Moreover (p' 1 ', cr*- 1 ^) solves @, twfeiie (p 1 - 2 ', ct*- 2 ^) satisfies (|3a[) /or all functions r\ £ C 1 (fi( 2 )) suc/i Wiat 
n = 0forx£ dD. (2) . 

Since the pair (P(i), o^i)) solves then p m is the mass density of an isolated single body supported 
in the region . Note however that we prove no results neither on the regularity of the boundary, nor 
on the regularity of the mass density in the interior of the body. These important open problems require 
further investigation. 

The identity (|28fl is an example of "Virial Theorem" , which in the context of particle mechanics gives 
information on how the energy of a system in equilibrium is distributed between kinetic and potential 
energy. Note finally that in view of (|15p , the energy of the minimizer tjj^ is negative, since ip(X) — X + d 
belongs to for 

Examples of stored energy functions 

The stored energy function of Ogden materials is given by 

L N 

w(X, F) = ^2ai(X) (Tr Cy- /2 + J2 b J ( x ) ( Tr Cof Cf 3 ^ + r ( det F) + h{X), (29) 

i=l 3=1 

where L, N £ N, C — F T ■ F is the (left) Cauchy-Green tensor, ji > 1, for i = 1, . . . L, Sj > 1, for 
j = 1, . . . N, r : (0, oo) ^> (0, oo) is a C 1 convex function such that 

F(z) > c\z~ a , for some ci > 0, 

di, bj : B — > (0, cxj) are measurable bounded functions such that 

d = essinfxss niin at{X) > 0, 

i—l,...L 

b = essinfxes niin bAX) > 0. 

i=l,...JV 

The function h is defined in such a way that the normalization condition 1)151) holds. As shown in |12l 
Th. 4.9-2], the stored energy function (|29|) is polyconvex and satisfies the coercive inequality in (w3)^' 
for 

p = max7j, q = maxj,-, (30) 

i j 

and a constant a > depending on d,b, ci. Moreover, assuming that 

|r'0)| < — (1+T(z)), for some c 2 >0, 
z 

the stored energy function (129[) satisfies (w4) as well, see [4] Sec. 2.4]. We remark that Ogden materials 
also meet other desirable physical requirements, such as material frame indifference, and provide a rich 
supply of case studies with important practical and theoretical applications [22] . 
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3 Proof of the main results 
3.1 Proof of Theorem CD 

To begin with we recall a result which will be used to show that the sets A^ 1 ' and A^ 2 ' are weakly closed. 
The proof can be found in [2], see also [121 Thm. 7.6-1]. 

Lemma 3. Let {tp n } C W 1 ' a (B) such that i}> n -± tp in W 1 ' a {B), Cof Vi/>„ H in L b (B) and det Vip H -± 
S in L C (B), for 

a > 2, a~ 1 +b~ 1 <l, c > 1. 

Then 

H = CofVtfj and 5 — det X7ip a.e. on B. 

In the following, the letter C will be used to denote various, possibly different, positive constants. 
Recall that @^{X,Y) = \tp(X) - ^(Y)^ 1 is a.e. defined on B x B, for all V £ A {i) , i = 1,2. 

Lemma 4. For all ip £ A^\ i = 1,2, and /or a!i < A < 3 ibc /lane 



0,/>(X, Y) A dX dY < C (^j det V 'ip (X)^ dX^j 



>B Jb 

Proof. By the the change of variables formula (|18[) . 

I^(B){x) ^(B){y) dxdy 



B JB 



Q^X^fdXdY = / / 



R 3j R 3 detWW- 1 (^))detV^(t/»- 1 (?/)) |a;-y| A ' 



By the Hardy-Littlewood-Sobolev inequality |18l Th. 4.3] 



9^(X, Y) A dXdY < C 



B JB 



U V(B) 



(detVVO °i>~ 1 
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L6-A , 



= C\ detV^ -1 ^))^ dx 

= c(J detViP(X)^e dX 

and the proof is complete. □ 
Using Lemma [4] with A = 1 and Holder's inequality we obtain 

i^m<o( f ^S^)' <c( f , ,3i, 



det Vip(X) 1 / 5 J ~ \J B detV^(X) 

where s is the exponent given in (w3) . Combining (|31f) with (w3)^' and the lower bound (|10p on p re f 
we have 

m > apo L detvmy - \M\vo»M*<m - c (J B det v*(xy ) " ■ ^ 

Since s > 1/3, we infer that I[ip] is bounded from below over A^\ The lower boundedness of I[ip] implies 
the existence of minimizing sequences in the space A^\ i.e. a sequence denoted by {ipn^} such that 

Ityn ] — ► m f ^> as n — ^ co. 

AM 

Owing to (|32[l. along any minimizing sequence we have 

[ d 4, < C, (33) 

Jb det V^\xy 
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and by pi[) we obtain that the sequence {E po t[ipn ']} is bounded. Since the sequence {ipn* 1 } is minimizing, 
then {E^ip^]} is also bounded. Using (w3) w , (jTOj and the bound detF < C|F| 3 , we infer that 

{V^i } is bounded in L P (B), 
{det Vi/'i } is bounded in L p/3 (B), 
{Cof Vip^} is bounded in L"(B). 

Finally, having prescribed the boundary value of {ipn^} C A^ and the average of {V>i 2) } C A {2 \ 
Poincare's inequality entails that 

{V-i } is bounded in W hP (B). 

It follows that any minimizing sequence possesses a subsequence {i/>n'} such that, for some ip^ G 
W^ P (B), S G L p/a {B), H G L q (B), 

mW 1 *(B), detV^^S mL p/ \B), CofV^^H in L q (B) 

and in addition, by the Rellich-Kondrasov imbedding theorem, 

V>i l) ^ (l \ in C°(E). 

Moreover, by Lemma [3] 

S = detVip (l \ H = Cof VV W - 

Now we show that 

V> (i) G.A (i) . (34) 

By weak convergence, det Vi/i' 1 ' > a.e. on 23. We claim that det W > a.e. on B. Indeed, suppose 
by contradiction that there exists a set L such that meas L > and det Vip^ = on L. Then, as shown 
in [121 pagg. 374-5], there exists a subsequence of {ipn* 1 } (not relabeled) which converges pointwise a.e. 
on L to a function ip^ satisfying det S/tp^ = on L. Consequently, by Fatou lemma, 

liminf _B s tr[^n ] = +oo, 

n— S-oo 

which implies that / = +00 for all ip £ However this is impossible, because the space is not 

empty, and so det W > holds a.e. on B. This completes the proof of (|34[1 when i = 1. Let us now 
consider the case i = 2. Clearly ^M 1 ' satisfies the last condition in the definition of A^ (by the weak 
convergence tp„ ip^ in L P (B)) and therefore it only remains to show that measS^(i) = 0. By (|18jl . 
the identity 

/ det VV4 15 (X)dX = meas (B) , 

JB 

holds for all minimizing sequences {ip 1 ^ 1 }. Passing to the limit n — > oo, using the weak convergence 
det Vtpn — det V^ 1 ' in L P ^ A {B) and the uniform convergence ipn — > ip^\ we obtain 



/ det X7ip (1) (X)dX = meas ip (1> (B) . 

JB 



Applying (|17[) to the left hand side we obtain 

card{(V' (1) )~ 1 (a:)}da; = meas ?/; (1) (B). (35) 



<il>W(B) 

If measS^fi) > 0, we would have 



/ card{(V' (1) )" 1 (x)}di = meas(V' <1) (S)\S'^(i))+ / card^ 1 ')^ (x)}dx 

> meas (ip^ (B) \S^(i)) + measS' l/ ,(i) = meastp^iB) 
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and therefore (|35[) leads to an absurd result. Hence meas S^cn = must hold, which completes the proof 
of igU for i = 2. 

Now, the assumptions (wl), (w2) imply 

liminf f Plsf (X)w(X,Vipn)dX> [ p lef (X)w(X, V^ (l) ) dX, 

n->oo J B J B 

see [2] and [121 Th. 7.7-1]. Thus the proof of existence of minimizers is complete if we show that 

£ P ot[Vi !) ] ->£ P ot[^ W ]. (36) 

We estimate 

\E pot [^]-E pot [^]\<C ( f \®,»{X,Y)-®, w {X,Y)\dXdY 



= C / 9 ii) e i , (i) \\rP^(X)-i>^(Y)\-\^(X)-^>(Y)\\dXdY 
Jbjb Vn 

<c f f e {i) e (i) \^(x)-^(x)\dXdY 

JBJB v ™ 

< cu^ - ^ 1^- ( B) ||e^w % w \w (BxB) 

< c'UV'i - v^'lli^cB)!!®^) llL2 (BxB )||e v ,(i)|| L 2 (BxB) . 

By Lemma[4] ({33j| , and since s > 1/2, {O^w} and {6^,<;)} are bounded in l?(B x 23); using the uniform 

convergence ip^ -> ip {i \ the claim J36J follows. 

To conclude our discussion on the existence of minimizers, we want to show that the convergence of 
the potential energy functional can also be obtained by using the properties of the spatial density along 
minimizing sequences, which is the argument used for fluids [24]. We define 

^4° = p^ )1 ^\by 

Before stating the next result, we remark that since {ipn^} is uniformly bounded in B, the supports of 
the spatial densities are all contained in a common compact region of R 3 . 

Proposition 1. For any minimizing sequence {ipn } C A^ z \ there exists a subsequence such that 

P+m "» L 1 ^), for all 1 < 7 < 1 + s. 

Moreover e po t[p^] -> e po t [£,/,(«)]■ 
Proof. First we observe that 



/ p (i) (x) 1+s dx=[ A ^*X' dX < \\p Tei H^ (B) 
Jr3 is det ViW, MX) s 



/s detV^ j (X) s ""~ w i B detVV>&°PO s 

1(M 3 ) 



whence by (|33[) . p (») is bounded in L 1+S (R 3 ). Since furthermore ||p ||l1(m 3 i = M, for all n, there 



exists G L 7 (R 3 ) and a subsequence {p^(i>} such that 



p^p^py inL^(R J ), K 7 <l + S . 
Moreover by the Dunford-Pettis Theorem, 

P^^pP inL^R 3 ). 
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Next we show that pi** 1 = p^(n almost everywhere. Let (f> £ C°(R' :i ). We have 



(pi l) - p 4 ,<i))<l>dx = lim / (p w - p^ w )<j>dx 



lim / p lei (X)[<f>(i>W{X))-<f>(i>W(X))}dX = 0, 



and the claim follows. Finally, the assertion e pot [p^(;)] — ► e po t[/f^,(o] follows by the well-known compact- 
ness properties of the functional e pot , see [201 pag. 125] and [231 Lemma 3.7]. 



3.2 Proof of Theorem H 

The claim £ L 1 (S7 < - 1 ^) follows by (w4) and the change of variable formula, for we have 

p le{ ((^)-\x)) 



f \a^(x)\dx = [ g(X,V^).(V^ W ) T 
JnW Jaw ot 



L 



X = (^( i >)- 1 (x) 



det V^WC^W)- 1 ^)) 



dx 



dw 
dF 



dX <K ( (l + w(X,Vip {i) ))dX < oo. 
Jb 



Next we show the validity of the identity (|28[) . Let r < 1 and set 

Clearly, i/^ 11 G yl (1) , for all r < 1. Moreover 

I[i>i 1) ]= f p le t(X)w(X, (1 - r)V^ (1) )dX + 1 £ P ot[^ (1) ] = Z(r). 
We claim that ?(•) is differentiable at r = 0. To prove this we estimate, for r > 0, 



w(X, Vip^) -w(X,Vip m )\ < 



j ITs fa*' (1 ~ sr ) v ^ (1) )) ds 



< r 



1 TiY||(X, (1 - S t)W (1) ) • (V^ (1) ) T ) 



< Kt(1+w(X, Vi/> (1) )) 
+ .JC / vv> t (1) ) - V^ (1) )|dt. 



Thus by Gronwall's inequality 

Vip^) - V^ (1) )| < Kre KT (l + w{X, Vip (1) )). 
The same estimate with r replaced by — r in the right hand side holds for r < 0, whence 
\w(X,Vipi 1) ) -w(X,Vip W )\ < K\r\e KlTl (1 + w(X,X7ip W )), Vr < 1. 

It follows that 

sup -L\w(X,S7i> { r 1) ) - w(X,V^ (1) )| < tfe^l+wpf,^ 1 ')). 
M<i/2 |r| 

By the Dominated Convergence Theorem, ?(•) is differentiable at r = and we have 
I'(0) = -f Prei{X) Tr V^ (1) ) • (W> (1) ) T ) + Bpot[^ (1) ]. 



□ 



(37) 
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Since tp^- 1 ' is a minimizer, then l'(0) = must hold: 

| B P rcf (X)Tr V^ (1) ) • (W> (1) ) T ) dX = £ pot [V» (1) ]. 

Upon the change of variable x = the latter transforms into (|28l) . 

Now let 77 (1) <E r> n( i)(R 3 ) and 

t? (2) : S^7 —>■ R 3 , t/ 2) G C^STt 3 )"), such that r? (2) (x) = 0, for x G <9S7 (2) . 
For t G R we define 

= </> W P0 +r7/(V Ci) (X)) + 2)-^ J PlotV (^(X))dX. 

We have 

V^ l) (X) = (J + rV I) ^ W (X))) ■ VV W (X), 

hence 

det = det V^ W det(J + rVr?^ (X))). 

Since det(I + tA) = 1 + rTr(A) + 0(t 2 ), there exists r > such that detVf/4^ > a.e. on B for 
|r| < r . Clearly ip ( r } £ A {2 \ for |r| < r . Moreover 

p rcf ^ 1) (X)dX = a. 

Thus £ .A (1) follows if we show that measS <i) = for |r| small enough. To this purpose we need 
the following result. 

Proposition 2. There exists C > such that 

\r) (i) (x) -rj (i) (y)\ < C\x-y\, for all x,y G , (38) 

i.e., r? w 6 C ' 1 ^). 

Proof. We shall use the following property of C 1 functions / : SI — s> R 3 defined on a regular domain SI: 
There exists a constant C, which depends only on fi, such that \f(x) — f(y)\ < C\\ V f\\i,°° (n) \x — y\, for 
all x, y G S7. When S7 is convex this is a trivial consequence of the Mean Value Theorem for vector valued 
functions (in this case we may choose C = 1); for general domains £7, the claim is proved as follows. 
A standard property of any domain S7 (see [121 pag. 224]) is that, for all x,y G fi, there exists a finite 
sequence of points Zi, ■ ■ ■ zjv+i, such that Z\ — x, zjv+i = y, Zi G Q, for alii = 2, . . .N, the open segment 
with end points Zi, Zi+i is all contained in S7, and the inequality 

JV 

^ |z* - 2j+i| < C|a; - y\ 
holds for a positive constant C which depends only on Q. Using this property we have 

N 

!/(*)- /(</)!= 



y ^[/(s2i + (l-s)zj+i)]ds 
< ||V/|| £ - (n) Y, - ^ C||V/|U- (n) |x - y|. 



Since Si' 2 ' is a regular domain, the inequality (|38|) for i = 2 is proved. For i = 1, let cf> £ Cc(R 3 ) such 
that 0| f (1) = j/ 1 ' and let 7? > such that S7*- 1 ' C -Br, where Br denotes the ball with center in the 
origin and radius R. Applying the preceding result we have 

\4>{x) - 4>{y)\ < C\x - y\, for all x, y G B R , 

where C = C(R). Restricting the previous inequality to x,y G SI*- 1 - 1 yields (|38[1 for i = 1. □ 
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Corollary 1. There exists r\ > such that Km) = K^ti), for \t\ < r\. In particular, meas-K^(i) = 
and thus 

meas S (i) = meas (7f (i> ) = 0, /or |r| < Ti. 

Proof. Clearly, -K^(i) C K <x>- Now, let A' 6 K^(X>- Thus there exists Y £ £> such that Y =fi X and 
4 1 \x) = ^ 1 \Y),i.e, 

^ W {X) + tt7(i/; (1) (X)) = ip {1) {Y) + t V (4> w {Y)). 

Using (|38p we obtain 

|i/> (1) (A) - > \4> (1) (X) - 001(1 - C*r). 

Thus ip (1) (X) = ip m {Y), for \r\ < (2C*)" 1 = n; hence X £ K^ m and K m C A^ (1) . □ 

In the following we shall assume that |r| < min{ro,ri}, so that ?/>)■ £ A' 1 '. Moreover, we simplify 
the notation by denoting both and r/ 2 ^ by rj. We continue by proving that the function r — > /[V>t '] 
is different iable at r = 0. We have 

lim / [^' ) ] ~ J ^ ( ' ) ] = i im i(E str [V>«] - E stI [^]) + lim i(£ pot [#] - E pot [^ (i) ]). (39) 

r— >0 7" t— >0 7" r— >0 T 

As shown by J. Ball in [3j pag. 13], the bound (w4) implies that the function r — I E B tr[ipT ] is differentiable 
at t — and 

lim-(E stI [^]^E stI [i,^])= C^-SstrCV^)) =/ (^-w(A,V^)) ProfdA 

= | e Tr ^(A, VV» W ) ■ (V^ W ) T ■ (Vr;(^W(A))) T ) p rcf dA. 
It will now be shown that a similar property holds for the potential term in (|39p . By (1381) . 

|^°(A) - V^'OO - W W (A) - V W 00)l < C\t\\^(X)-^(Y)\; 

thus the estimate 

|^(A)-^«(y)| > |V (!) (A)-^>(F)|/2, i.e., e^)(A,V) <29^,(A,y), 
holds for \r\ small enough. Moreover 

|9^)(A,r) - %^(X,Y)\ = Q^ ) (X,Y)e^ ) (X,Y)\\^\X)~^\Y)\ - \^{X) - ^{Y)\\ 

< e^ ) (X,Y)Q^ i) (X,Y)\(\^(X)-^(X)\ + \^)(Y) - ^\Y)\) 

<c|r|e 4i) (A,y)e^ w (A,y). 

Thus for \r\ small enough we have 



By LemmalU and since s > 1/2, the function in the right hand side is integrable on B x B. Then, by the 
Dominated Convergence Theorem, and since 

lim 1 r« (XV\ A fyy^ (^W(A) - ^W(Y)) ■ (^ w (*)) - ^OO)) 
for almost all X,Y € B, we obtain 



lim iCEpottV^] " £ P ot[V (!) ]) = -lim^ /" /" p rof (A)p rcf (F) (e^ (A, y) - 9^, (A, y)) 



dXdY 



p rcf (A) Prof (y) ^) {X )-^)(yw 
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Substituting into (|39|) and using that f// 1 ' is a minimizer yields 

+ Is Is p - l(x)p ^ Y) WH^WWW = 



Upon the change of variable x = ip^(X), y = V OOj the latter equation transforms into (|3aJ, and the 
proof of Theorem [5] is concluded. 
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